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The production and the observability of A(1520), K°(892) 
$ and A(1232) hadron resonances in central Pb+Pb collisions 
at 160 AGeV is addressed. The rescattering probabilities of 
the resonance decay products in the evolution are studied. 
Strong changes in the reconstructable particle yields and spec- 
tra between chemical and thermal freeze-out are estimated. 
Abundances and spectra of reconstructable resonances are 
predicted. 



Strange particle yields and spectra are key probes to 
study excited nuclear matter and to detect the transi- 
tion of (confined) hadronic matter to quark-gluon-matter 
(QGP) 10-0- The relative enhancement of strange and 
multi-strange hadrons, as well as hadron ratios in central 
heavy ion collisions with respect to peripheral or proton 
induced interactions have been suggested as a signature 
for the transient existence of a QGP-phase 0. 

Unfortunately, the emerging final state particles re- 
member relatively little about their primordial source, 
since they had been subject to many rescatterings in 
the hadronic gas stage This has given rise to 

the interpretation of hadron production in terms of ther- 
mal/statistical models. Two different kinds of freeze-outs 
are assumed in these approaches: 

1. a chemical freeze-out, were the inelastic flavor 
changing collisions processes cease, roughly at an 
energy per particle of 1 GeV [03 , 

2. followed by a later kinetic/thermal freeze-out were 
also clastic processes have come to an end and the 
system decouples. 

Chemical and thermal freeze-out happen sequentially at 
different temperatures (T ch w 160 - 170 MeV, T th « 
120 MeV) and thus at different times. 

To investigate the sequential freeze-out in heavy ion 
reactions at SPS the Ultra-relativistic Quantum Molec- 
ular Dynamics model (UrQMD 1.2) is applied 
UrQMD is a microscopic transport approach based on 
the covariant propagation of constituent quarks and di- 
quarks accompanied by mesonic and baryonic degrees of 
freedom. It simulates multiple interactions of ingoing 
and newly produced particles, the excitation and frag- 
mentation of color strings and the formation and decay 
of hadronic resonances. At present energies, the treat- 
ment of sub- hadronic degrees of freedom is of major im- 
portance. In the UrQMD model, these degrees of free- 



dom enter via the introduction of a formation time for 
hadrons produced in the fragmentation of strings [p^|-p^[ . 
The leading hadrons of the fragmenting strings contain 
the valence-quarks of the original excited hadron. In 
UrQMD they are allowed to interact even during their 
formation time, with a reduced cross section defined by 
the additive quark model, thus accounting for the original 
valence quarks contained in that hadron jl3|,[l4]]. Those 
leading hadrons therefore represent a simplified picture of 
the leading (di)quarks of the fragmenting string. Newly 
produced (di)quarks do, in the present model, not in- 
teract until they have coalesced into hadrons - however, 
they contribute to the energy density of the system. For 
further details about the UrQMD model, the reader is 
referred to Ref. HIO). 

Let us start by asking whether a microscopic non- 
equilibrium model can support the ideas of sequential 
chemical and thermal break-up of the hot and dense mat- 
ter? To analyze the different stages of a heavy ion colli- 
sion, Fig. [j] shows the time evolution of the elastic and in- 
elastic collision rates in Pb+Pb at 160AGeV. The inelas- 
tic collision rate (full line) is defined as the number of col- 
lisions with flavor changing processes (e.g. 7T7t — > KK). 
The elastic collision rate (dashed dotted line) consists of 
two components, true elastic processes (e.g. 7T7t — > nn) 
and pseudo-elastic processes (e.g. ttt: — > p — > 7T7t). While 
elastic collision do not change flavor, pseudo-elastic col- 
lisions arc different. Here, the ingoing hadrons are de- 
stroyed and a resonance is formed. If this resonance de- 
cays later into the same flavors as its parent hadrons, this 
scattering is pseudo-elastic. 

Indeed the main features revealed by the present mi- 
croscopic study do not contradict the idea of a chemical 
and thermal break-up of the source as shown in Fig. |l| 
(top) However, the detailed freeze-out dynamics is much 
richer and by far more complicated as expected in sim- 
plified models: 

1. In the early non-equilibrium stage of the AA colli- 
sion (t < 2 fm/c) the collision rates are huge and 
vary strongly with time. 

2. The intermediate stage (2 fm/c< t < 6 fm/c) is 
dominated by inelastic, flavor and chemistry chang- 
ing processes until chemical freeze-out. 

3. This regime is followed by a phase of dominance of 
elastic and pseudo-elastic collisions (6 fm/c< t < 
11 fm/c). Here only the momenta of the hadrons 
change, but the chemistry of the system is mainly 



1 



unaltered, leading to the thermal freeze-out of the 
system. 

Finally the system breaks- up (t > 11 fm/c) and the scat- 
tering rates drop exponentially. 

Fig. |] (bottom) depicts the average energj|j] per par- 
ticle at midrapidity (\y — y cm \ < 0.1). One clearly ob- 
serves a correlation between chemical break-up in terms 
of inelastic scattering rates and the rapid decrease in en- 
ergy per particle. Thus, the suggested phenomenological 
chemical freeze-out condition of ss 1 GeV/particle is also 
found in the present microscopic model analysis. 




t cms (fm/c) 



FIG. 1. Topdnelastic and (pseudo-)elastic collision rates in 
Pb+Pb at 160AGeV. r c h and 7th denote the chemical and 
thermal/kinetic freeze-out as given by the microscopic reac- 
tion dynamics of UrQMD. Bottom: Average energy per par- 
ticle at midrapidity (\y — y C m\ < 0.1) as a function of time. 

To verify this scenario, we exploit the spectra and 
abundances of A(1520), if (892) and other resonances 



To compare to thermal model estimates, the energies are 
calculated from interacted hadrons with the assumption p 2 z = 
(pl+p 2 y )/2. This assures independence of the longitudinal 
motion in the system and the chosen rapidity cut. 



which unravel the break-up dynamics of the source be- 
tween chemical and thermal freeze-out. In the statisti- 
cal model interpretation of heavy ion reactions the reso- 
nances are produced at chemical freeze-out. If chemical 
and thermal freeze-out are not separated - e.g. due to an 
explosive break-up of the source - all initially produced 
resonances are reconstructable by an invariant mass anal- 
ysis in the final state. However, if there is a separation 
between the different freeze-outs, a part of the resonance 
daughters rescatter, making this resonance unobservable 
in the final state. Thus, the relative suppression of res- 
onances in the final state compared to those expected 
from thermal estimates provides a chronometer for the 
time period between the different reaction stages. Even 
more interesting, inelastic scatterings of the resonance 
daughters (e.g. Kp — > A) might change the chemical 
composition of the system after 'chemical' freeze-out by 
as much as 10% for all hyperon species. While this is 
not in line with the thermal/statistical model interpreta- 
tion, it supports the more complex freeze-out dynamics 
encountered in the present model. 

To answer these questions we address the experimen- 
tally accessible hadron resonances: At this time <£> and 
A(1520) have been observed in heavy ion reactions at 
SPS energies 00 following a suggestion that such a 
measurement was possible pC] . SPS p9[ and RHIC ex- 
periments [J2ll report measurement of the if (892) sig- 
nal, _and RHIC has already measured both the K and 
the K°. In the SPS case, the A(1520) abundance yield is 
about 2.5 times smaller than expectations based on the 
yield extrapolated from nucleon-nucleon reactions. This 
is of highest interest in view of the A enhancement by fac- 
tor 2.5 of in the same reaction compared to elementary 
collisions. 

As an explanation for this effective suppression by a 
factor 5, we show that the decay products (7r, A, etc.) 
produced at rather high chemical freeze-out temperatures 
and densities have rescattered. Thus, their momenta do 
not allow to reconstruct these states in an invariant mass 
analysis. However, even the question of the existence of 
such resonance states in the hot and dense environment 
is still not unambiguously answered. Since hyperon res- 
onances are expected to dissolve at high energy densities 
(see .e.g. p2fl) it is of utmost importance to study the 
cross section of hyperon resonances as a thermometer of 
the collision. 

The present exploration considers the resonances, 
A(1232), A*(1520), K°(892) and $. The properties of 
these hadrons are depicted in Tabic |. 



Particle 


Mass (MeV) 


Width (MeV) 


r (fm/c) 


A(1232) 


1232 


120 


1.6 


A* (1520) 


1520 


16 


12 


K°(892) 


893 


50 


3.9 


$(1020) 


1019 


4.43 


44.5 



TABLE I. Properties of investigated resonances. 
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Fig. U shows the rapidity densities for A(1232), 
A*(1520), K°(892) and $ in Pb(160 .4GeV)Pb, b< 3.4 fm 
collisions. Fig. |^ (left), shows the total amount of decay- 
ing resonances. Here, subsequent collisions of the decay 
products have not been taken into account - i.e. when- 
ever a resonance decays during the systems evolution it 
is counted. However, the additional interaction of the 
daughter hadrons disturbs the signal of the resonance 
in the invariant mass spectra. This lowers the observable 
yield of resonances drastically as compared to the primor- 
dial yields at chemical freeze-out. Fig. ||(right) addresses 
this in the rapidity distribution of those resonances whose 
decay products do not suffer subsequent collisions - these 
resonances are in principle reconstructable from their de- 
cay products. Note that reconstructable in this context 
still assumes reconstruction of all decay channels, includ- 
ing many body decays. 
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FIG. 2. Rapidity densities for A(1232), A* (1520), K°(892) 
and $ in Pb(160 ,4GeV)Pb, b< 3.4 fm collisions. Left: All 
resonances as they decay. Right: Reconstructable resonances. 



At SPS energies, the rapidity distributions dN/dy may 
by described by a Gaussian curve 



d^ (y)=Axex H"^ 



(i) 



with parameters given in Table || However, in this ap- 
proximation the details of the rapidity distributions are 
lost. Especially the dip in the rapidity distributions of 
the A(1520) is not accounted for. 

The rescattering probability of the resonance decay 
products depends on the cross section of the decay prod- 
uct with the surrounding matter, on the lifetime of the 



surrounding hot and dense matter, on the lifetime of the 
resonance and on the specific properties of the daugh- 
ter hadrons in the resonance decay channels. This leads 
to different 'observabilities' of the different resonances: 
Rescattering influences the observable $ and A* yields 
only by a factor of two, due to the long life time of those 
particles. In contrast strong effects are observed in the 
K* and A yields, which are suppressed by more than a 
factor three. 

One can use the estimates done by |g3j in a statistical 
model, and try to relate the result of the present micro- 
scopic transport calculation to thermal freeze-out param- 
eters. The ratios of the 4ir numbers of reconstructable 
is A(1520)/A = 0.024 and K°(892)/K+ = 0.25. In terms 
of the analysis by (23) , the microscopic source has a life- 
time below 1 fm/c and a freeze-out temperature below 
100 MeV. Thus, the values obtained from UrQMD seem 
to favor a scenario of a sudden break-up of the initial 
hadron source, in contrast to the time evolution of the 
chemical and thermal decoupling as shown in Fig. [l|. 
However, note that these numbers are based on the above 
mentioned thermal scenario. In fact, fitting hadron ratios 
of UrQMD calculations for central Pb+Pb interactions 
at 160 AGeV with a Grand Canonical ensemble yields a 
chemical freeze-out temperature of 150 — 160 MeV p4| . 

In fact, the rescattering strength depends on the phase 
space region studied. Fig. || addresses the longitudinal 
momentum distribution of the rescattering strength. The 
probability R to observe the resonance H is given by 



R 



H — ► hih 2 (reconstructable) 
H fti/i 2 (all) 



(2) 
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FIG. 3. Rapidity dependent ratio R of reconstructable res- 
onances over all resonances of a given type as they decay. For 
A(1232), A* (1520), K°(892) and $ in Pb(160 AGeV)+Pb, 
b< 3.4 fm collisions. 
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The 'observability' decreases strongly towards central 
rapidities. This is due to the higher hadron density at 
central rapidities which increases the absorption proba- 
bility of daughter hadrons drastically. Unfortunately, in 
the case of the meson this decrease of the observable 
yield increases the discrepancies between data and mi- 
croscopic model predictions (e.g. Q). 

The absorption probability of daughter hadrons is not 
only rapidity dependent but also transverse momentum 
dependent. The decay products are rescattered preferen- 
tially at low transverse momenta. Fig. |] depicts the in- 
variant transverse momentum spectra of reconstructable 
A(1232), A*(1520), K°(892) and $ in Pb(160 AGeV)Pb, 
b< 3.4 fm collisions. 



Fig. [5] directly addresses the p t dependence of the ob- 
servability of resonances. The present model study sup- 
ports a strong p t dependence of the rescattering proba- 
bility. This effects leads to a larger apparent tempera- 
ture (larger inverse slope parameter) for resonances re- 
constructed from strongly interacting particle. A similar 
behavior has been found for the $ meson in an indepen- 
dent study by [ p5| . This effective heating of the <E> spec- 
trum might explain the different inverse slope parameter 
measured by NA49 ((f) — > K + K~) as compared to the 
NA50 (<j) -> n+tL-) collaboration 
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FIG. 4. Transverse momentum spectra at \y — y cm | < 1 of 
reconstructable resonances for A (1232), A* (1520), K°(892) 
and $ in Pb(160 ^GeV)Pb, b< 3.4 fm collisions. The lines 
are to guide the eye. 



Particle (All) 


A 


a 


A 


181 


1.49 


A* 


0.89 


1.20 


K* 


22.0 


1.23 


$ 


2.22 


1.09 


Particle (Rec.) 


A 


a 


A 


46.3 


1.60 


A* 


0.42 


1.27 


K* 


6.95 


1.34 


$ 


1.62 


1.10 



TABLE II. Gaussian fit parameters for the rapidity densi- 
ties of all and reconstructable resonances. 
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FIG. 5. Ratio R of reconstructable resonances over all reso- 
nances of a given type as a function of transverse momentum. 
For A(1232), A* (1520), K°(892) and $ in Pb(160 AGeV)Ph, 
b< 3.4 fm collisions. 

In conclusion, central Pb+Pb interactions at 160 
AGeV are studied within a microscopic non-equilibrium 
approach. The calculated scattering rates exhibit signs 
of a chemical and a subsequent thermal freeze-out. The 
time difference between both freeze-outs is explored with 
hadronic resonances. The rapidity and transverse mo- 
mentum distributions of strange and non-strange res- 
onances are predicted. The observability of unstable 
(strange) particles (A(1520), <fi, etc.) in the invariant 
mass analysis of strongly interacting decay products is 
distorted due to rescattering of decay products from 
chemical to thermal freeze-out. Approximately 25% of 
</>'s and 50% of A* 's are not directly detectable by recon- 
struction of the invariant mass spectrum. The rescatter- 
ing strength is strongly rapidity dependent. Rescattering 
of the decay products alters the transverse momentum 
spectra of reconstructed resonances. This leads to higher 
apparent temperatures for resonances. Inelastic collisions 
of anti-Kaons from decayed strange resonance alter the 
chemical composition of strange baryons by up to 10%. 
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events 1400 

A/event 654.46 

A* /event 2.60 

K° (892) /event 66.2 

$/event 5.917 



TABLE III. 4-7T yields of all decaying resonances in 
Pb(160AGeV)+Pb, b< 3.4 fm. 



y dn/dy(A) dn/dy(A(1520)) dn/dy(K°(892)) dn/dy($) 

0.125 169.12 0.78 21.36 2.16 

0.375 167.40 0.89 20.36 2.05 

0.625 163.89 0.79 19.20 1.89 

0.875 155.86 0.7 17.56 1.65 

1.125 143.33 0.58 14.99 1.33 

1.375 128.78 0.52 12.83 1.10 

1.625 110.98 0.36 10.31 0.76 

1.875 90.976 0.26 7.554 0.48 

2.125 67.511 0.16 4.815 0.26 

2.375 47.786 0.09 2.378 0.09 

2.625 36.04 0.03 0.802 0.03 

2.875 21.583 0.01 0.2 0.00 

3.125 5.3613 0.00 0.028 

3.375 0.2719 0.00 

3.625 0.0017 

3.875 

TABLE IV. Rapidity density of all decaying resonances in 
Pb(160AGeV)+Pb, b< 3.4 fm. 



events 1400 

A/event 178.9 

A*/event 1.28 

K° (892) /event 22.61 

$/event 4.35 



TABLE V. 47T yields of all reconstructable resonances in 
Pb(160AGeV)+Pb, b< 3.4 fm. 



y 


dn/dy(A) dn/dy(A(1520)) dn/dy(K°(892)) dn/dy(<I>) 


0.125 


44.30 


0.348 


6.44 


1.601 


0.375 


43.03 


0.437 


6.59 


1.484 


0.625 


42.46 


0.37 


6.14 


1.342 


0.875 


40.05 


0.332 


5.75 


1.21 


1.125 


37.08 


0.292 


5.16 


1.008 


1.375 


33.74 


0.262 


4.47 


0.797 


1.625 


29.82 


0.200 


3.81 


0.579 


1.875 


25.72 


0.154 


2.91 


0.361 


2.125 


21.14 


0.087 


2.14 


0.22 


2.375 


16.57 


0.050 


1.18 


0.0752 


2.625 


12.43 


0.018 


0.46 


0.022 


2.875 


8.570 


0.007 


0.12 


0.002 


3.125 


2.755 


0.001 


0.02 





3.375 


0.165 





0.00 





3.625 


0.001 











3.875 















TABLE VI. Rapidity density of reconstructable resonances 
in Pb(160AGeV)+Pb, b< 3.4 fm. 



0.125 


169.12 


0.78 


21.36 


0.375 


167.40 


0.89 


20.36 


0.625 


163.89 


0.79 


19.20 


0.875 


155.86 


0.7 


17.56 


1.125 


143.33 


0.58 


14.99 


1.375 


128.78 


0.52 


12.83 


1.625 


110.98 


0.36 


10.31 


1.875 


90.976 


0.26 


7.554 


2.125 


67.511 


0.16 


4.815 


2.375 


47.786 


0.09 


2.378 


2.625 


36.04 


0.03 


0.802 


2.875 


21.583 


0.01 


0.2 


3.125 


5.3613 


0.00 


0.028 


3.375 


0.2719 





0.00 


3.625 


0.0017 








3.875 
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